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^ Abstract 

' Context. This paper reports the results obtained on the photometric redshifts measurement and accuracy, and cluster tomography in 

Q \ the ESO Distant Cluster Survey (EDisCS) fields. 

Aims. We present the methods used to determine photometric redshifts to discriminate between member and non-member galaxies 
and reduce the contamination by faint stars in subsequent spectroscopic studies. 

Methods. Photometric redshifts were computed using two independent codes both based on standard spectral energy distribution 
(SED) fitting methods (Hyperz and G. Rudnick's code). Simulations were used to determine the redshift regions for which a reliable 
determination of photometric redshifts was expected. The accuracy of the photometric redshifts was assessed by comparing our 
estimates with the spectroscopic redshifts of ~ 1400 galaxies in the 0.3 < z < 1.0 domain. The accuracy expected for galaxies fainter 
than the spectroscopic control sample was estimated using a degraded version of the photometric catalog for the spectroscopic sample. 
Results. The accuracy of photometric redshifts is typically <x( Az/(1 + z)) ~ 0.05 ± 0.01, depending on the field, the filter set, and 
the spectral type of the galaxies. The quality of the photometric redshifts degrades by a factor of two in cr(Az/(l + z)) between the 
brightest (/ 1 22) and the faintest (/ ~24-24.5) galaxies in the EDisCS sample. The photometric determination of cluster redshifts in 
the EDisCS fields using a simple algorithm based on z pAo , is in excellent agreement with the spectroscopic values, such that Sz ~0.03- 
0.04 in the high-z sample and Sz ~0.05 in the low-z sample, i.e. the z p i w t cluster redshifts are at least a factor ~ (1 + z) more accurate 
than the measurements of z p i m , for individual galaxies. We also developed a method that uses both photometric redshift codes jointly 
to reject interlopers at magnitudes fainter than the spectroscopic limit. When applied to the spectroscopic sample, this method rejects 
~ 50 - 90% of all spectroscopically confirmed non-members, while retaining > 90% of all confirmed members. 
Conclusions. Photometric redshifts are found to be particularly useful for the identification and study of clusters of galaxies in large 
surveys. They enable efficient and complete pre-selection of cluster members for spectroscopy, allow accurate determinations of the 
cluster redshifts based on photometry alone, and provide a means of determining cluster membership, especially for bright sources. 

Key words. Galaxies: clusters — Galaxies: distances and redshifts — Galaxies: photometry — Galaxies: evolution 
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1. Introduction 

Photometric redshifts are becoming an important tool in cosmo- 
logical studies based on large and/or deep photometric surveys. 
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Different studies have been devoted to the detailed analysis of 
photometric redshift accuracy in different contexts (e.g. Ilbert et 
al. 2006, Feldmann et al. 2006, Mobasher et al. 2007, Banerji et 
al. 2008, Margoniner & Wittman 2008, Hildebrandt et al. 2008, 
Ilbert et al. 2009). The robust evaluation of the accuracy reached 
by photometric redshifts requires homogeneous deep photomet- 
ric data and a large dataset of spectroscopic redshifts for the 
same field. Simulations can be used to achieve uniform coverage 
in parameter spaces beyond the limits of spectroscopic surveys, 
in particular when missing information about certain redshift do- 
mains and/or spectroscopic types of galaxies. 

Several papers used a recalibration between data and tem- 
plate models to improve the precision of photometric redshifs 
(e.g. Coe et al. 2006, Ilbert et al. 2006, Feldmann et al. 2006, 
Mobasher et al. 2007, Capak et al. 2007, Ilbert et al. 2009), a 
method that requires a large and representative training set of 
spectroscopic redshifts. However, model templates, optimized to 
achieve the highest possible accuracy in a given catalog/field, are 
not necesarily optimal in all cases because systematic problems 
in the catalog photometry could remain unrecognized during the 
template calibration process. A more robust estimate of photo- 
metric redshifts accuracy can be achieved for large datasets ac- 
quired in different independent fields. This is the approach used 
in this paper. 

The ESO Distant Cluster Survey (hereafter EDisCS) is an 
ESO Large Programme designed to study the evolution of cluster 
galaxies over a significant fraction of cosmic time (White et al. 
2005). The 20 clusters included in the EDisCS sample were se- 
lected from the Las Campanas Distant Cluster Survey (LCDCS; 
Gonzalez et al. 2001) with redshifts ranging between ~0.4 and 1. 
More details about the survey and the cluster selection procedure 
can be found in the paper by White et al. (2005), and the EDisCS 
websitf]. The EDisCS programme includes homogeneous and 
deep photometry with ESO VLT and NTT (optical and near-IR; 
White et al. 2005, Aragon-Salamanca et al. in preparation) and 
multi-object spectroscopy with ESO VLT (Halliday et al. 2004, 
Milvang-Jensen et al. 2008), as well as other follow-up observa- 
tions with HST/ACS (Desai et al. 2007), narrowband Ha imag- 
ing (Finn et al. 2005) and XMM data (Johnson et al. 2006). 

This paper is also intended to be the reference for the pho- 
tometric redshifts and the cluster membership criteria adopted 
by the EDisCS collaboration, and used in the different EDisCS 
papers dealing with cluster membership and related quantities 
(e.g. De Lucia et al. 2004, White et al. 2005, Clowe et al. 2006, 
Poggianti et al. 2006, De Lucia et al. 2007, Desai et al. 2007, 
Rudnick et al. 2009). Photometric redshifts are particularly use- 
ful when used with cluster/structure finding algorithms, because 
they help to ensure that time-consuming spectroscopic observa- 
tions are optimized (e.g. Li & Yee 2008). 

The paper is organized as follows. In Sect. [2l we summarize 
the characteristics of the relevant photometric and spectroscopic 
data. A technical description of the photometric redshift methods 
and related procedures is provided in Sect. [3] The photometric 
redshift accuracy in the EDisCS survey is addressed in Sect. [4] 
using three different approaches: 1) simulations are used to de- 
termine the redshift regions for which a reliable determination of 
photometric redshifts is expected; 2) the actual quality achieved 
in this survey is estimated by direct comparison between the 
photometric and spectroscopic redshifts, and 3) the accuracy ex- 
pected for galaxies fainter than the spectroscopic control sample 
is estimated using a degraded version of the spectroscopic sam- 
ple catalog. Section [5] presents the comparison between spec- 
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troscopic and photometric determinations of cluster redshifts, as 
well as the results obtained on cluster tomography in the EDisCS 
fields. The photometric cluster membership criteria adopted by 
the EDisCS collaboration is introduced and discussed in Sect. [6] 
Discussion and conclusions are given in Sect. [7] The follow- 
ing cosmological parameters are adopted throughout this paper: 
Q.A = 0.7, Q m = 0.3, and H = 70 km s' 1 Mpc~ l . Magnitudes 
are given in the Vega system. 

2. Photometric and spectroscopic data 

We use the ground-based photometric catalogs and spectro- 
scopic redshifts obtained by the EDisCS collaboration for 20 
clusters of galaxies with redshifts ranging between 0.4 and 1.0 
(White et al. 2005). Although the final redshift distribution of 
this sample is found to be fairly uniformly distributed within this 
redshift interval, the original filter set was designed to bracket 
the relevant wavelength domain at z~0.5 (low-z sample) and 
z~0.8 (high-z sample). Photometric redshifts and related quan- 
tities depend strongly on the wavelength domain covered by the 
photometric Spectral Energy Distributions (hereafter SED), i.e. 
the filter set. Throughout the paper, we therefore retain the orig- 
inal division of the clusters into the "low-z" and the "high-z" 
samples. 

Deep optical photometric data was acquired with FORS2 at 
the VLT, in B VI and VRI bands for the low-z and the high-z clus- 
ter samples, respectively. The photometric depth (5cr in 1 "radius 
aperture) is typically 26.4(B), 26.2(V) and 24.8(1) in the low-z 
sample, and 26.5(V), 26.0(R) and 25.2(1) in the high-z sample 
(see also White et al. 2005, Table 1). The field of view covered 
by these data is 6.5' x 6.5'. Seeing conditions were excellent 
during all imaging observations, ranging typically between 0.5" 
and 0.8" (see White et al. 2005 for details). Deep near-IR im- 
ages were also obtained for almost all clusters with SOFI at 
the NTT, in K, and JK S for the low-z and the high-z samples, 
respectively (details are provided by Aragon-Salamanca et al. 
2009, in preparation). The photometric depth (5<x) is typically 
22.8 in J and 21.5 in K s . These data cover a field of 4.2'x 6.0' 
at low-z, and 4.2'x 5.4' at high-z. Photometry was performed 
on seeing-matched images using SExtractor v.2.2.2 (Bertin & 
Arnouts 1996). Table[3]summarizes the filter set used to compute 
photometric redshifts for each cluster in the EDisCS sample. 

Spectroscopic data in the EDisCS fields were obtained 
during three observing runs with FORS2 at VLT, using the 
600RJ+19 grism. The wavelength domain covered by our ob- 
servations ranged between ~ 5300 and 9000 A. More details can 
be found in the reference papers by Halliday et al. (2004) and 
Milvang-Jensen et al. (2008). The total number of good qual- 
ity spectra acquired per field ranged between ~60 and 100 for 
the low-z sample, and was around ~100 for the high-z sample. 
There are typically 30-50 confirmed members in every cluster. 
We use objects with either secure spectroscopic redshifts (here- 
after type 1) or medium quality, slightly tentative redshifts (here- 
after type 2) to characterize the behavior of photometric redshifts 
and cluster-membership criteria. Objects with tentative redshifts 
(~ 50% secure, type 3) represent less than 2% of the total sam- 
ple and are mostly used for illustration pusposes. The total num- 
ber of spectroscopic redshifts available is 637(977) in the low- 
z(high-z) samples, from which the total number of secure (secure 
+ slightly tentative, i.e. type 1+2) redshifts in the 0.3 < z < 1.0 
domain are 544(564) and 854(885) respectively for the low-z 
and high-z samples (see also TableHJi. 

An important issue when deriving photometric redshifts for a 
given galaxy is to construct its SED using magnitudes and corre- 
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sponding fluxes derived for identical aperture sizes in each of the 
bandpass images. Photometric SEDs were obtained from seeing- 
matched images according to the following scheme. For "iso- 
lated" objects (SExtractor flag =0), photometric redshifts were 
derived from isophotal magnitudes measured within the refer- 
ence I-band isophotal region corresponding to 1.5 cr of the local 
background noise. In the case of "crowded" objects (SExtractor 
flag >0), we used magnitudes computed within 1" radius aper- 
tures. This scheme enabled us to improve the SED determination 
in crowded regions, while increasing the S/N for isolated galax- 
ies. 

We did not use the standard SExtractor errors because these 
are known to underestimate the error for dithered data where ad- 
jacent pixels are correlated. We determined our errors instead 
by means of a set of empty aperture simulations as described in 
White et al.(2005). The accurate determination of the errors is 
important because photometric redshifts are sensitive to photo- 
metric errors. 

A correction for Galactic extinction was also included for 
each cluster field according to the E(B-V) derived from Schlegel 
et al. (1998) for the center of the cluster. The E(B-V) correc- 
tions in the EDisCS fields typically ranged between 0.03 and 
0.08 magnitudes (see TableO. 

Bright unsaturated stars were used as secondary standards 
to check the consistency of our photometric system for deriv- 
ing photometric redshifts and, when required, to introduce small 
zero-point corrections in the photometric catalogs. In practice, 
we compared the color-color diagrams of observed stars in our 
fields with the expected positions derived using the Pickles li- 
brary (Pickles 1998). Stars at this stage are objects selected with 
SExtractor stellarity index > 0.95 that belong to the stellar se- 
quence in the I-band isophotal-radius versus aperture-magnitude 
diagrams. This procedure was particularly useful during the first 
photometric runs to correct near-IR imaging data for the ef- 
fects of non-photometric conditions. Due to the addition of high- 
quality observations, the zero-point corrections improved suc- 
cessively during the lifetime of the project, in addition to the 
quality of the photometric catalogs and related quantities, such 
as photometric redshifts. The results presented here were ob- 
tained with the final version of the photometric catalogs, for 
which the zero-point corrections are negligible, apart from two 
cases: C11138-1133 (AV =0.10 and AJ =-0.15), and C11232- 
1250 (AJ =-0.20). 

All the results published by the EDisCS collaboration since 
2004 were obtained with the current and final version of the pho- 
tometric catalogs used in this paper, publically available from 
the EDisCS websitfl The last version of EDisCS photometric 
redshifts was obtained in April 2006. The quality of this final 
version with respect to the previous ones (since 2004) is about 
the same in terms of accuracy (i.e. systematic offsets, dispersion 
and catastrophic failures; see criteria in Sect. 0). The main dif- 
ferences come from the related quantities which are provided 
in addition to photometric redshifts (e.g. absolute magnitudes, 
photometric classification of galaxies, ...). 



3. Photometric redshifts 

Photometric redshifts (hereafter z p hot ) were computed using two 
different codes: a modified version of the public code Hyper^ 
(Bolzonella et al. 2000), and the code of Rudnick et al. (2001), 
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with the modifications introduced by Rudnick et al. (2003) (here- 
after GR code). The two codes use different approaches based on 
SED fitting procedures, as summarized below. The reader is re- 
ferred to the reference papers for a more detailed description of 
the codes themselves. Here we summarize only the main rele- 
vant settings and modifications. 

Hyperz results were initially used by the EDisCS collabo- 
ration for three main purposes: to determine a first guess for 
each cluster redshift, to help in spectroscopic pre-selection, and 
to reduce the contamination by faint stars during spectroscopic 
observations. Subsequently, the two codes were jointly used 
to establish cluster membership in magnitude-limited samples 
using their respective normalized probability distributions (see 
Sect. I6). 

We used 14 galaxy templates with Hyperz: 

- Eight evolutionary synthetic SEDs computed with the 2003 
version of the Bruzual & Chariot code (Bruzual & Chariot 
1993, 2003), spanning a grid of ages between 0.0001 and 
13.5 Gyrs, with Chabrier (2003) IMF and solar metallicity 
(a delta burst -SSP-, a constant star-forming system, and 6 
r-models with exponentially decaying SFR). 

- A set of 4 empirical SEDs compiled by Coleman, Wu and 
Weedman (1980) (hereafter CWW) to represent the local 
population of galaxies, with fixed age, extended to wave- 
lengths A < 1400 A and A > 10000 A using the equivalent 
Bruzual & Chariot spectra. 

- Two starburst galaxies (SB1 and SB2) from the Kinney et al. 
(1996) template library. 

The internal reddening law is taken from Calzetti (2000), and 
considered as a free parameter with Ay ranging between and 
1.5 magnitudes (E(B-V) between and ~ 0.45 mags). When an 
object is not detected in a given filter, the flux in this filter is set 
to zero with an error bar corresponding to the limiting magnitude 
that corresponds to a S/N ratio ~ 1 in this filter. Absolute magni- 
tudes Mb for galaxies were allowed between -24 < Mb < —11, 
which imposed a relatively weak prior on z p hat while prevent- 
ing obvious catastrophic identifications in the case of degenerate 
solutions. P(z) were normalized within the permitted redshift in- 
terval according to this restriction. 

The GR code is based on the non-negative linear combina- 
tion of redshifted galaxy templates, for which: 

- The set of 4 CWW empirical templates described above. 

- The starburst galaxies SB1 and SB2 from the Kinney et al. 
(1996). 

- A 10 Myr old, single stellar population burst obtained from 
the 1999 version of the Bruzual & Chariot (1993) code with 
Salpeter (1955) IMF and solar metallicity. 

The 4 CWW and starburst templates were extended from 
their published short-wavelength limits (1400 and 1232 A re- 
spectively) by a power-law fit to the 1400-1800 and 1240-1740A 
wavelength ranges, respectively, using Bruzual & Chariot (2003) 
spectra. 

There are no limitations on absolute magnitude in the GR 
code, and the direct flux measurements were used for all galax- 
ies, even when an object was not formally detected. 

Hyperz z p hot were computed in the range < z < 6, whereas 
GR ones span the < z < 2 range. The upper limit had a 
negligible impact on the z p hot value itself and related quantities, 
except in the normalization of the z p hot probability distribution 
(P(z)). When deriving the cluster membership criteria presented 
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in Sect. [6] P(z) computed with the two codes were normalized 
within the < z < 2 interval, and Hyperz Zphot were also re- 
stricted to this interval. 

3. 1 . Photometric discrimination between galaxies, stars and 
quasars 

The method used to discriminate photometrically between 
galaxies, stars, and quasars, was based on Hyperz, and closely 
followed the developments presented in Hatziminaoglou et al. 
(2000). For stars, it was based on a standard SED fitting mini- 
mization about z ~ using the complete library of stellar tem- 
plates by Pickles (1998). The galactic E(B-V) correction was 
considered as a free parameter, ranging between and the cor- 
responding value for the field given in Table [3] For quasars, we 
used a library of synthetic spectra similar to Hatziminaouglou et 
al. (2000), and the same prescriptions as for galaxies, apart from 
the absolute magnitude limitation. 

In practice, the usual Hyperz z p hot for galaxies and quasars, 
and the best fit with the stellar library were computed for each 
object, and three classification parameters were given to quan- 
tify the goodness of the best fit as a galaxy/star/quasar (re- 
spectively Nc, N, and Nq). The object was "rejected" as a 
galaxy/star/quasar when its x 1 excluded it at higher than the 
95% confidence level (N=0). The object was "fully compatible" 
when the probability associated with the reduced x 1 exceeded 
90% (N=2). The object was "undetermined" (N= 1) in all the 
other intermediate cases. This classification allowed us to de- 
fine different samples of objects in these fields, either galaxies 
(with Nc > 1, irrespective of the star type) or stars (with N* > 
1 and Nc<l). These classification criteria were used during the 
spectroscopic runs to lower the contamination due to faint and 
red stars to values below ~ 10% (Halliday et al. 2004; Milvang- 
Jensen et al. 2008), and they are also used below in § 15.31 The 
quasar classification was not considered for the spectroscopic 
preselection. Johnson et al. (2006) used this classification to 
identify possible AGNs detected in X-rays in these fields. 

3.2. Photometric classification of galaxies 

Galaxies were classified into five different spectral types, cor- 
responding to their rest-frame photometric SED: (1) E/S0, (2) 
Sbc, (3) Scd, (4) Im, and (5) SB (starbursts). These types corre- 
spond to the simplest empirical templates given above, namely 
the four SEDs compiled by CWW for the local galaxies, plus the 
Kinney starbursts. This classification corresponds to the best fit 
templates for the GR code. In the case of Hyperz, it is the best fit 
of the rest-frame SED at z p hot • The classification obtained with 
the two codes is in perfect agreement, excepted for catastrophic 
identifications. We use this classification in § H to address the 
Zphot accuracy as a function of the spectral type. 

3.3. Photometric redshift catalogs 

Photometric redshifts for EDisCS catalogs computed with the 
two different codes are publically available from the EDisCS 
websit^]. 

EDisCS online catalogs also include an optimized flag for 
the discrimination between galaxies and stars, based on the com- 
bination between the above Hyperz criterion, a similar fit using 
the GR code (flagcfl=l when the object is well fit as a star) (1), 
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z(model) 

Figure 1. Photometric versus model redshifts retrieved from 
BVIK, SEDs (low-z sample), for 10 5 simulated galaxies uni- 
formly distributed between < z < 2. The diagram displays the 
number density of galaxies in linear scale. 

the SExtractor stellarity flag (flag,s£ A ) (2), and the size for bright 
objects (3), i.e. 

1. {N, >1 ANDN G <1 }ORflag OT = 1 

2. flag S£ , > 0.95 

3. r/,(I) < T thresh if l tot <22.5 

where r/,(I) is the half-light radius in the I band, and r t kresh 
is the threshold radius determined from the stellar locus in the 
corresponding iy, versus total I magnitude. 

4. Determination of the photometric redshift 
accuracy 

The expected accuracy of z P hot as a function of redshift depends 
strongly on the photometric accuracy and the filter set used to 
derive the photometric SEDs. As explained in §[2] and in White 
et al. (2005), we introduced an accurate determination of pho- 
tometric errors to address the former issue. The filter set used 
by EDisCS contains a relatively small number of filters because 
it was designed originally to cover the relevant wavelength do- 
main in the rest frame of low z~0.5 and high z~0.8 clusters. 
In particular, they bracketed the 4000 A break for the relevant 
redshift intervals, but they were not designed to explore the full 
redshift domain. For this reason, we address the z p hoi accuracy in 
three different ways described below. Using simulations, we first 
determine the redshift ranges for which we expect a reliable esti- 
mate of Zphot for the low and high-z filter sets, and the ideal (max- 
imum) accuracy expected from simple SED fitting. Secondly, the 
quality of EDisCS Zphot achieved is estimated by direct compar- 
ison with the spectroscopic redshifts. In a third subsection, we 
determine the z p hot accuracy expected for galaxies fainter than 
the spectroscopic control sample. 
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Figure 2. Photometric versus model redshifts retrieved from 
VRIJK, SEDs (high-z sample), for 10 5 simulated galaxies uni- 
formly distributed between < z < 2. The diagram displays the 
number density of galaxies in linear scale. 



The following quantities were computed to quantify the 
Zphot accuracy, where z S pec stands for both "spectroscopic" and 
"model" redshifts: 

- The systematic deviation between z p hot and z spe c > (A z ) = 
2 A z /N, given by the mean difference between these two 
quantities, where A z = z spe c - Zphot and N is the total num- 
ber of galaxies. 

- The standard deviation cr z = VIX4: ~ ( A z)) 2 /( N ~ 1)> ex- 
cluding catastrophic identifications, defined here in a con- 
servative way for those galaxies with |A Z |= \z sp ec - Zphot I > 

0. 3X(1+ Z sp ec )• 

- The median absolute deviation cr zMAD = 1 .48 x median 
I z sp ec - Zphot \, which is less sensitive to outliers. 

- The normalized median absolute deviation defined as 
cr(Az/(l + z)) = 1.48 x medianQzspec - z p hot \/(l+z sp ec )), 

- The percentage of catastrophic identifications (/%), i.e. 
galaxies "lost" from their original spectroscopic redshift bin, 
with|A,| >0.3x(l+z !FC ). 

- The percentage of galaxies included in a given photomet- 
ric redshift interval that are catastrophic identifications (g%), 

1. e. galaxies that contaminate the sample because they are in- 
correctly assigned to the redshift bin. 

4. 1 . Expected accuracy from simulations 

Photometric redshift determinations are based on the detection 
of strong spectral features, such as the 4000 A break, Lyman 
break or strong emission lines. In general, broad-band filters al- 
low only detection of strong breaks and are insensitive to the 
presence of emission lines, apart from when the contribution of 
a line to the total flux in a given filter is higher than or simi- 
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Figure 3. These figures present the expected z P hot accuracy de- 
rived from simulations, for galaxies with photometric quality 
similar to the spectroscopic sample. The top panel displays the 
systematic deviation (A z ) ± cr z , the standard deviation excluding 
catastrophic identifications, for the BVIK, (low-z, black dots) 
and VRIJK, (high-z, open/red dots) filter sets, within the relevant 
redshift domain. The botton panel displays a plot of cr(Az/(l +z)) 
as a function of redshift for the same filter combinations. 



lar to the photometric errors, as happens in the case of AGNs 
(Hatziminaoglou et al. 2000). 

To determine the redshift domains where a reliable measure- 
ment of Zphot can be obtained given the filter sets used in the low 
and high-z samples, we completed a series of simulations assum- 
ing a homogeneous redshift distribution in the redshift interval 
< z < 2. These simulations were performed using Hyperz re- 
lated software, and z p hot computed with both Hyperz and the GR 
code, but the results should be representative of the general be- 
havior of pure SED-fitting z p hot codes. Synthetic catalogs contain 
10 5 galaxies within this redshift interval, for each filter set, span- 
ning all the basic spectral types defined in § 13.21 with uniform 
redshift distribution. Photometric errors in the different filters 
were assigned following a Gaussian distribution with <x scaled 
to magnitudes according to Am 2.51og[l + l/(S/N)], where 
S IN is the signal-to-noise ratio, which is a function of apparent 
magnitude. Here we used the mean S /N achieved in the different 
filters for the spectroscopic (bright) sample, i.e. for objects rang- 
ing between 1=18.5(19.0) and 22.0 in the low-z (high-z) sample, 
and the same settings used for z p hot computation on real catalogs. 
In this way, the results obtained from simulations should be con- 
sidered to be "ideal" but still consistent with those derived in 
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Table 1. z p hot accuracy from simulations for galaxies in the spectroscopic sample, based on Hyperz. 



Clusters 


Redshift interval 


Galaxy type 


<A Z > 




°~ z.MAD 


or(Az/(l +z)) 


1% 


g% 


Low-z 


0.3 < z phB , < 1 


all 


0.061 


0.131 


0.124 


0.076 


2.2 


6.0 


High-z 


0.3 < Zphot < 1 


all 


0.015 


0.101 


0.074 


0.045 


6.6 


2.8 


High-z 


0.45 < z phot < 1 


all 


0.019 


0.108 


0.080 


0.046 


8.2 


2.8 


Low-z 


0.3 < Zphot < 1 


E/S0 


-0.036 


0.121 


0.087 


0.052 


3.4 


2.2 




Sbc 


0.075 


0.108 


0.176 


0.105 


2.2 


5.0 






Scd 


0.110 


0.082 


0.176 


0.107 


0.8 


1.9 






Im 


0.023 


0.089 


0.086 


0.053 


0.4 


5.6 






SB 


0.095 


0.124 


0.126 


0.075 


0.6 


11.6 


High-z 


0.3 < Zphot < 1 


E/S0 


-0.019 


0.040 


0.031 


0.019 


0.0 


0.2 




Sbc 


0.037 


0.081 


0.077 


0.048 


0.1 


0.7 






Scd 


0.030 


0.072 


0.070 


0.043 


0.0 


1.4 






Im 


0.050 


0.077 


0.085 


0.050 


0.2 


1.1 






SB 


-0.022 


0.127 


0.087 


0.054 


7.4 


3.5 



Notes - The accuracy reached for the different spectral types of galaxies is also presented. The information given is the same as in TableH for the 
current spectroscopic sample. 

Table 2. z p hot accuracy expected for the faintest galaxies in the EDisCS sample, in the low-z and high-z fields, based on Hyperz. 



Clusters 


Redshift interval 


Galaxy type 


<A ; > 




C z.MAD 


o-(Az/(l +z)) 


1% 


g% 


Low-z 


0.3 < z p i, ot < 1 


all 


-0.003 


0.153 


0.158 


0.105 


4.0 


1.0 


High-z 


0.3 < zphot < 1 


all 


-0.037 


0.159 


0.156 


0.095 


2.2 


6.2 


High-z 


0.45 < z p hot < 1 


all 


-0.027 


0.156 


0.151 


0.091 


1.0 


6.9 


Low-z 


0.3 < z ph „, < 1 


E/S0 


-0.025 


0.155 


0.172 


0.113 


3.2 


0.0 




Sbc 


0.025 


0.151 


0.158 


0.106 


4.6 


0.7 






Scd 


0.021 


0.129 


0.124 


0.082 


2.7 


1.8 






Im 


-0.012 


0.172 


0.161 


0.113 


8.8 


2.5 






SB 


-0.096 


0.156 


0.238 


0.155 


6.7 


5.6 


High-z 


0.3 < Zphot < 1 


E/S0 


-0.040 


0.155 


0.151 


0.091 


2.0 


2.6 




Sbc 


-0.041 


0.164 


0.172 


0.108 


2.1 


1.9 






Scd 


-0.017 


0.166 


0.146 


0.092 


1.7 


8.7 






Im 


-0.026 


0.143 


0.142 


0.088 


0.9 


2.8 






SB 


-0.075 


0.150 


0.191 


0.125 


6.3 


23.7 



Notes - The accuracy reached for the different spectral types of galaxies is also presented. The information given is the same as in TablefU for the 
current spectroscopic sample. 



§ 14.21 Because of the limited scope of these simulations, we do 
not intend to explore all possible domains of parameter space, 
but focus instead on studying the main systematics introduced 
by the photometric system. 

Figures Q]and E] display the photometric versus model red- 
shifts retrieved from BVIK, (low-z sample) and VRIJK, (high-z 
sample), respectively. Figure [3]and Table [TJ summarize the qual- 
ity of Zphot in the simulations using Hyperz, within the relevant 
redshift domain. The same results are presented in Table [9] for 
the GR code. 

Some systematic trends are clearly visible in figures [TJ |2]and 
[3] The high-z filter combination provides a higher quality and 
smaller systematic deviations than the low-z one. This trend was 
expected because of the more complete and contiguous spectral 
coverage of the high-z set. The lack of an R filter for the low-z 
sample introduces a systematic trend in A- ~ 0.05-0.08 at z > 0.3; 
the highest Zphot quality for this sample is expected to be around 
cr(Az/(l + z));0.06 at 0.4 <z< 0.6, i.e. within the sensitive red- 
shift domain. Because of the lack of B-band photometry for the 
high-z sample, the highest quality results are expected at z > 0.4. 
This trend is indeed observed in the simulations. However, the 
quality achived for high-z simulated data at z < 0.4 with Hyperz 
is expected to be overestimated with respect to real data, because 
templates and models are drawn from the same parent set. This 
is a general criticism of simulations used in assessing the real- 



istic performance of z p i w t quality. Also z p hot quality depends on 
the spectral type of galaxies. With respect to the average qual- 
ity presented in Fig. [3]for a uniform distribution of types, early 
types exhibit up to a ~50% improvement in cr(Az/(l + z)) with 
Hyperz in the redshift domains where the filter sets bracket the 
4000 A break. 

The results for the GR code in Table [9] are similar in average 
to Hyperz 's ones. The quality tend to be slightly better for the 
bluest spectral types, whereas it is worse for early types. This 
trend can be explained by the broad paramer space spanned by 
the simulations, the same used by Hyperz for z p i w t determina- 
tions, as compared to the GR code. The noise for late-type galax- 
ies in Hyperz tend to be dominated by degeneracies, whereas 
the GR code cannot properly fit highly-reddened E-Sbc galaxies. 
Note that the uniform distribution in redshift and types in these 
simulations does not represent a realistic population of galaxies. 

4.2. Comparison with spectroscopic redshifts 

The Zpiwt accuracy achieved for EDisCS was estimated by di- 
rect comparison with its 1449 spectroscopic redshifts in the 
0.3 < z < 1 interval (type 1+2). Although spectroscopic tar- 
gets in our sample were strongly biased in favor of cluster mem- 
bers with Zphot within the interval z c i usrer + 0.2 (see Halliday et 
al. 2004; Milvang-Jensen et al. 2008, and § I5H , the geometri- 
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Figure 4. (z spe c - z p hot ) versus spectroscopic redshift for the low-z clusters in the EDisCS sample. Solid (red) circles, open (blue) 
circles, and crosses correspond to objects with good (type 1), medium (type 2) and tentative (type 3) spectroscopic redshift deter- 
minations, respectively. Error bars in z P hot correspond to lcr. Dot-dashed lines display Zspec - Zphot = ±0.1 to guide the eye. (see text 
for details). 



cal configuration of slits in one hand, and the need for reference 
field galaxies on the other hand, ensured that non-member galax- 
ies were also targeted during the spectroscopic runs. In principle, 
these field galaxies should allow us to extend the present study to 
the < Zspec < 1 interval, according to the restrictions imposed 
by the set of filters (see § 14. 11 1. 



Table [3] presents the z P hot accuracy obtained with Hyperz 
and G. Rudnick's (GR) code for the different low-z and high- 
z clusters, based on the direct comparison with the spectro- 
scopic sample, excluding stars. Only type 1 objects, i.e. objects 
with secure spectroscopic redshifts, were taken into account. In 
Table [3] systematic deviations and <x, were computed over all 
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z(spectroscopic) z(spectroscopic) 

Figure 5. (z spec - z p hot ) versus spectroscopic redshift for the high-z clusters in the EDisCS sample. Solid (red) circles, open 
(blue) circles, and crosses correspond to objects with good (type 1), medium (type 2) and tentative (type 3) spectroscopic redshift 
determinations, respectively. Error bars in z p hot correspond to lcr Dot-dashed lines display z spe c - z p hot = ±0. 1 to guide the eye. (see 
text for details). 



the 0.3 < Zspec ^ 1 interval. In the low-z bin, the two clusters 
with only B VI photometry were excluded from the sample when 
computing the average z P kot quality over the cluster sample. In 
the high-z bin, we excluded Cll 122-1 136 from the cluster statis- 
tics, because of the small number of z sp ec available in this field. 

Our main result is that there is no significant systematic shift, 
neither in the (A z ) ± cr z nor in the cr(Az/(l + z)) results, with 



respect to the values expected from simulations in § 14.11 with 
some field-to-field differences discussed below. The accuracy of 
Zphot ranges usually between cr(Az/(l + z)) ~ 0.05 for Hyperz 
and ~ 0.06 for GR, both for the low-z and the high-z samples. 
This result is in good agreement with the highest possible ac- 
curacy expected from ideal simulations in the low-z case, it is ~ 
25% worse than ideal expectations in the high-z case for Hyperz, 
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and compatible with GR results for late type galaxies (early-type 
errors were overestimated, as mentioned in Sect. 14. It . The dif- 
ferential trend between low and high-z samples with respect to 
simulations in both codes can be explained because of the dif- 
ferent population of "bright" galaxies in these fields, low-z and 
high-z samples containing a smaller and larger fraction of late- 
type galaxies respectively, compared to the uniform average pop- 
ulation in simulated data (e.g. De Lucia et al. 2007). This effect 
is clearly seen in Table |4] The fraction of objects lost from (or 
spuriously assigned to) the relevant redshift interval according 
to the definitions given above (/% and g%), is negligible in the 
low-z sample and typically below 5% for the high-z one. 

Table |4] summarizes the results obtained for the entire low-z 
and high-z samples, for both type 1 (secure) and type 1 + type 2 
(both secure and slightly tentative redshifts) spectroscopic data. 
The results are similar in both cases. In the high-z bin, z p hot ac- 
curacy improves slightly when the sample is restricted to the 
0.45 < Zspec ^ 1 interval, where both Hyperz and GR codes yield 
the same cr(Az/(l + z)) ~ 0.052. This effect is easily understood 
because at z;0.45 the rest-frame 4000 A break is found short- 
ward of the V-band filter. Table [4] also summarizes the z p hot ac- 
curacy achieved for the different spectral types of galaxies in 
the entire sample, i.e. all type 1, 2 and 3 spectra. Tentative type 
3 galaxies represent less than 2% of the total sample, and the 
results remain unchanged whith respect to type 1+2. The num- 
ber of galaxies as a function of the spectral type given in this 
table correspond to Hyperz. Excepted for catastrophic identifi- 
cations, the classification obtained with the two codes is in per- 
fect agreement. The lowest quality cr(Az/(l +z)) values are mea- 
sured as expected for the bluest galaxy types (SB), for both the 
Hyperz and GR codes. Early types display the highest quality 
results with Hyperz, whereas GR code has lower quality results 
for early types in the high-z sample. This trend may indicate that 
the CWW templates provide an inappropriate description of the 
SEDs of early types at intermediate redshifts. 

The comparison between the Hyperz and GR codes, either on 
a cluster-by-cluster basis or as a function of the filter combina- 
tion, yields similar results, even though these codes are based on 
different approaches and have different strengths/weaknesses. In 
general, Hyperz results are found to be of slightly higher accu- 
racy than GR's ones (by < 20% in <x(Az/(l + z))), but both are in 
close agreement with the expectations under "ideal" conditions. 
An interesting trend is that the quality of both codes is highly 
correlated, in the sense that the highest and the lowest quality 
results (in terms of cr( Az/(1 + z)) and systematics) are found for 
the same clusters. Given the homogeneous photometry of the 
EDisCS project, this trend can hardly be explained by the use of 
an incomplete or imperfect template set (as suggested by Ilbert 
et al. 2006), because in such a case we would expect the same 
systematic behavior in all fields, given a certain filter set, as is 
observed in § 14. II with Hyperz. In contrast, different systematic 
trends are observed for different clusters, which are then found 
to be almost identical for the two independent Zphot codes. This 
behavior suggests strongly that the origin of the systematics is 
more likely to be the input photometric catalog rather than the 
Zphot codes and templates. In particular, we cannot exclude small 
remaining zero-point shifts in our data, approximately equal to 
or less than ~ 0.05 magnitudes, because we are limited by the 
accuracy of the stellar templates (see Sect.|2]i. 

A brief discussion of particular aspects of z p hot accuracy in 
the low-z and high-z samples is given below. 



4.2.1. Low-z cluster fields 

Figure [4] displays a direct comparison between the spectro- 
scopic and the photometric redshifts for the low-z clusters in the 
EDisCS sample. Hyperz was used to derive z p hot in this figure, 
but the results with the GR code are very similar, as discussed 
above. Error bars in z P hot correspond to a lcr confidence level in 
the photometric redshift probability distribution P(z), i.e. to the 
68% confidence level computed through the Ay 2 increment for a 
single parameter (Avni 1976). Figure [TBI shows the comparison 
between spectroscopic and photometric redshifts for the entire 
low-z sample, obtained with Hyperz and GR codes, as well as 
the z S pec - Zphot distribution. 

The Zphot quality in this sample ranges usually between 
0.04 < cr(Az/(l + z)) < 0.07 with both Hyperz and GR codes, with 
some exceptions. On the one hand, Cll 119-1 129 and C11238- 
1 144 were only observed in BVI, which produces lower quality 
Zphot and a higher fraction of catastrophic identifications. These 
two clusters were not included when deriving the mean values 
in Table [3] On the other hand, CI 1232- 1250 was observed in J 
in addition to BVIK, and this provides a more accurate z P hot es- 
timate with respect to average with Hyperz, although there is no 
clear improvement with the GR code. Compared to simulations, 
the systematic trend A, ~ 0.05-0.08 at z > 0.3 is far smaller in real 
data, whereas <x(Az/(l + z)) is in agreement with ideal results. 

4.2.2. High-z cluster fields 

Figure [5] displays a direct comparison between the spectro- 
scopic and the photometric redshifts for the high-z clusters in 
the EDisCS sample. Error bars in z p hot correspond to lcr confi- 
dence level in the photometric redshift probability distribution 
P(z). Figure [T7] shows the comparison between spectroscopic 
and photometric redshifts for the whole high-z sample, obtained 
with Hyperz and GR codes, as well as the z spec -Zphot distribution. 

The Zphot quality in this sample usually ranges between 
0.05 < cr(Az/(l + z)) < 0.08 with both Hyperz and GR codes, with 
some exceptions. The statistics in Cll 122-1 136 is based on a 
small number of spectroscopic redshifts, hence we exclude this 
cluster when deriving the mean values in Table[4] Two out of the 
ten clusters in the high-z sample are actually in a redshift range 
typical of the low-z sample. Indeed, in the case of C11037-1243a 
and Cll 138-1 133, the low redshift of the cluster implies that B- 
band photometry is required to ensure that an accurate z P hot mea- 
surement is achieved, although the quality of their z P hot measure- 
ments is close to average. As seen in Fig. [5] individual error bars 
are larger in these two fields than in the other high-z clusters. 
Compared to simulations, there is no systematic trend in A z as 
expected, whereas cr(Az/(l + z)) is in agreement with ideal re- 
sults. 



able 3. Summary of z P hot accuracy achieved for the individual low-z and high-z clusters. 





Hyperz 


GR 


cluster (Low-z) 


Filters 


^cluster 


E(B-V) 


total 


NC1 ^ 


/A S 




& z.MAD 


(T{l\Z/\i + ZD 


1% 


g% 


/A \ 




°~ zMAD 


u ynz/ yi < 


z>) 


1% 


g /o 


(1) 




(3) 


(4) 


(5) 


(6) 


17) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 




(17) 
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H 1 01 8-1 71 1 


BVIK 
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0773 


63 


56 
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1 8 





cll059-1253 


BVIK, 


0.456 


0.0332 


78 


67 


0.018 


0.064 


0.069 


0.045 


0.0 


2.9 


0.015 


0.082 


0.069 


0.047 




0.0 


1.5 


CU119-1129 


BVI 


0.550 


0.0332 


57 


45 


0.061 


0.117 


0.133 


0.090 


11.1 


0.0 


0.032 


0.089 


0.103 


0.066 




11.1 


0.0 


cl 1202- 1224 


BVIK, 


0.424 


0.0583 


59 


46 


-0.007 


0.111 


0.113 


0.079 


0.0 


0.0 


0.041 
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0.0 


0.0 


C11232-1250 


BVIJK, 


0.541 


0.0596 


93 


77 


-0.013 


0.083 


0.060 


0.039 


1.3 


2.7 


0.022 
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0.097 


0.066 




0.0 


2.5 
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11 


8 
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0.072 
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0.072 


0.0 


0.0 
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12.5 


0.0 
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78 
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0.0 
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0.0 
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62 


55 
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0.0 
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68 


63 
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0.0 


0.0 


0.048 
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0.089 


0.058 




0.0 


0.0 
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BVIK,, 


0.496 


0.0821 


68 


58 


0.007 


0.088 


0.069 


0.046 


0.0 


0.0 


-0.001 


0.093 


0.091 
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0.0 


0.0 
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C11040-1155 


VRIJK, 


0.704 


0.0505 


116 


104 


0.028 


0.068 


0.080 


0.046 


0.0 


1.9 


0.021 


0.072 


0.076 


0.046 




0.0 


1.9 


C11054-1146 


VRIJK, 


0.697 


0.0374 


102 


94 


0.025 


0.090 


0.090 


0.054 


0.0 


3.3 


0.020 


0.127 


0.114 


0.068 




2.1 


4.3 


cl 1054- 1245 


VRIIK,, 


0.750 


0.0378 


98 


77 


-0.001 


0.087 


0.094 


0.055 


0.0 


5.8 


0.063 


0.071 


0.111 


0.064 




0.0 


7.0 


clll03-1245b 


VRIJK, 


0.703 


0.0481 


98 


86 


0.012 


0.088 


0.085 


0.052 


2.3 


4.1 


-0.004 


0.121 


0.112 


0.069 




7.0 


8.4 


CU122-1136 


VRIJK, 


0.640 


0.0404 


12 


7 


-0.012 


0.105 


0.180 


0.107 


14.3 


14.3 


-0.005 


0.115 


0.061 


0.042 




0.0 


0.0 


ell 138-1133 


VRIJK, 


0.479 


0.0274 


110 


94 


-0.002 


0.088 


0.088 


0.056 


0.0 


1.1 


-0.035 


0.080 


0.080 


0.054 




0.0 


2.0 


C11216-1201 


VRIJK, 


0.794 


0.0449 


116 


96 


0.034 


0.079 


0.096 


0.053 


1.0 


7.8 


0.040 


0.096 


0.136 


0.077 




0.0 


7.6 


C11227-1138 


VRIJK, 


0.635 


0.0468 


110 
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-0.030 


0.115 


0.087 


0.054 


1.0 


2.0 


-0.004 
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0.075 


0.047 




0.0 


3.0 
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0.762 


0.0793 


104 


92 


-0.018 


0.108 


0.090 


0.053 


3.3 


2.2 


0.010 


0.107 


0.086 


0.052 




1.1 


1.1 


Average 












0.003 


0.093 


0.091 


0.055 


0.8 


3.2 


0.011 


0.096 


0.099 


0.060 




1.3 


4.1 


rms 












0.021 


0.015 


0.007 


0.007 


1.1 


2.2 


0.027 


0.019 


0.020 


0.010 




2.1 


2.6 



otes - The information given in this table is: (1) cluster identification, (2) filter set used to compute z,,hot , (3) spectroscopic z chlsKr , (4) E(B-V) at the cluster center, (5) total number of redshifts 
mailable in the field, (6) total number of spectroscopic redshifts used for accuracy determinations, and z p hot accuracy achieved with Hyperz and GR codes: (7)(13) systematic deviation between 
,kot an d z S pec , (8)(14) standard deviation o%, (9)(15) median absolute deviation {(T zMAD ), (10)(16) normalized median absolute deviation (<x(Az/(l + z))), (11)(17) fraction of catastrophic identi- 
;ations (/%), (12)(18) fraction of spurious identifications (g%). Clusters excluded in the computation of the averaged accuracy are given in italic (see text). The spectroscopic z c i uster in column 
.) corresponds to the most prominent cluster in the field. Other clusters have been identified in cll037 (z =0.578), clll03 (z =0.959 and 0.626), clll38 (z =0.455), cll227 (z =0.583), cll301 
=0.397), and cll354 (z =0.595) (see Milvang- Jensen et al. 2008 for details). 



able 4. Summary of z P hot accuracy achieved for the whole low-z and high-z samples with Hyperz and GR codes. 
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0.054 


0.9 


3.3 


0.012 


0.100 


0.094 


0.058 


1.3 


4.0 




1 + 2 


all 


885 


0.004 


0.100 


0.090 


0.055 


1.1 


3.3 


0.010 


0.104 


0.095 


0.058 


1.4 


3.9 


High-z 


< 7 , , < l 
u.tj ^ iphot — 1 




all 


736 


0.013 


0.087 


0.087 


0.051 


0.7 


3.6 


0.024 


0.086 


0.088 


0.053 


0.4 


4.1 




1 +2 


all 


754 


0.014 


0.088 


0.087 


0.052 


0.9 


3.5 


0.023 


0.089 


0.089 


0.053 


0.5 


4.0 


Low-z 


0.3 < < 1 


1 + 2 + 3 


E/SO 


232 


0.006 


0.071 


0.069 


0.045 


0.0 


0.0 


0.036 


0.085 


0.091 


0.061 


0.0 


0.0 




1 + 2 + 3 


Sbc 


167 


0.027 


0.128 


0.104 


0.069 


1.2 


1.3 


0.036 


0.104 


0.092 


0.066 


0.6 


1.3 






1 + 2 + 3 


Scd 


113 


0.027 


0.103 


0.075 


0.047 


2.7 


0.0 


-0.004 


0.102 


0.077 


0.047 


4.4 


0.9 






1 + 2 + 3 


to 


46 


-0.007 


0.104 


0.087 


0.056 


().() 


4.0 


0.051 


0.166 


0.102 


0.062 


0.0 


0.0 






1 + 2 + 3 


SB 


18 


0.018 


0.153 


0.159 


0.094 


11.1 


5.0 


-0.019 


0.204 


0.167 


0.105 


5.6 


5.2 


High-z 


0.3 < z p h,t < 1 


1 + 2 + 3 


E/SO 


263 


0.040 


0.070 


0.088 


0.052 


0.8 


0.8 


0.025 


0.114 


0.126 


0.077 


0.8 


1.9 




1 + 2 + 3 


Sbc 


256 


-0.009 


0.103 


0.089 


0.055 


1.6 


1.6 


0.005 


0.108 


0.110 


0.068 


1.6 


2.7 






1 + 2 + 3 


Scd 


181 


0.017 


0.094 


0.094 


0.056 


1.1 


1.7 


-0.010 


0.079 


0.063 


0.038 


2.8 


7.8 






1 + 2 + 3 


to 


112 


-0.012 


0.079 


0.075 


0.047 


0.0 


0.9 


0.024 


0.066 


0.075 


0.044 


0.0 


0.0 






1 + 2 + 3 


SB 


81 


-0.083 


0.135 


0.173 


0.113 


2.5 


19.4 


0.012 


0.137 


0.107 


0.070 


2.5 


10.0 



otes - Results for type 1 and type 1+2 spectroscopic data are presented in this table, as well as the accuracy reached for the different spectral types of galaxies in the whole sample (i.e. all type 1, o 

and 3 redshifts). The information given is the same as in Table [3] I 
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4.3. Expected accuracy for galaxies fainter than the 
spectroscopic sample 

We determine the z P hot accuracy expected for galaxies fainter 
than the spectroscopic control sample used in § 14.21 i.e. galax- 
ies with magnitudes typically ranging between 1=18.5(19.0) and 
22.0 in the low-z (high-z) sample, in particular for the faintest 
galaxies in the EDisCS sample. The concept is to derive z P hot on 
a degraded version of the photometric catalog for the spectro- 
scopic sample (in terms of S /N), using the same recipes and 
settings as the main catalogs. This method was preferred instead 
of simulations because it uses the observed SEDs of the control 
sample instead of an arbitrary mixture of spectral types at a given 
redshift. 

Degradated catalogs were generated from the original (spec- 
troscopic) ones, to reproduce the photometric properties of the 
faintest galaxies in the EDisCS sample. The mean / magnitude 
was set to be < / >=24.00(24.5) for the low-z(high-z) clus- 
ter fields, corresponding to a S /N ~ 5. For all the other j fil- 
ters, magnitudes were scaled according to the original SEDs, i.e. 
keeping colors unchanged: m new {j) = m(j) + [24.0/24.5 - /] 
Photometric errors as a function of apparent magnitudes were 
introduced and assigned as in § 14.11 In this case, Aml ew (j) = 
[2.5 log[l + 1 /(S /N)]] 2 -A 2 m(j), where Am(j) is the catalog error 
corresponding to m(j), and S /N = S/AT(mQ))10- a4(m ^>-' nW) . 
This procedure conserves globally the colors of galaxies. The 
main caveat is the fact that this noisy population does not nec- 
essarily match the true color distribution of the faintest galaxies 
in the sample. However, it is useful to estimate the degradation 
expected in z p kot accuracy between the brightest and the faintest 
galaxies because of the lowered S/N. Because z P hot quality is 
quite insensitive to spectroscopic quality, we added the type 1 + 
type 2 spectroscopic catalogs. 

Table [2] and [10] summarize the results obtained for the 
faintest galaxies using Hyperz and the GR code respectively. 
These tables can be compared directly with Table [4] The qual- 
ity of photometric redshifts degrades typically by a factor of two 
in <x( Az/(1 + z)) between the brightest (7;22) and the faintest 
(/ ~24-24.5) galaxies in the EDisCS sample. Most of the trends 
observed in Table|4]for the spectroscopic sample are found in the 
Table [2] for the simulations of the faintest sample, in particular 
the lack of systematics in A z , and the higher quality results in the 
high-z bin. The fraction of catastrophic identifications increases, 
but remains typically below ~ 5%. The difference in z p hot qual- 
ity between early and late types is smaller for the faintest galaxy 
sample. In this case, the simulation results with the GR code are 
found to be of slightly higher accuracy than Hyperz 's ones (by 
<2Q% in <x(Az/(l +z))). In Sect. 16.31 we comment on the impli- 
cations that this results will have for the calculation of member- 
ship using the photometric redshifts. 

5. Photometric determination of cluster redshifts 

5.1. Spectroscopic sample preselection 

Before the first spectroscopic runs, cluster redshifts were esti- 
mated from the first z P kot catalogs using Hyperz. Spectroscopic 
targets were selected mainly to have z p hot within the interval 
^duster ± 0.2, or absolute P(z c i llster ) > 0.5, and according to the 
magnitude selection (see Halliday et al. 2004). Although the re- 
sults discussed in this Sect, were obtained with Hyperz, they 
should be representative of the general behavior of all SED- 
fitting Zphot codes. 

The photometic cluster redshifts were computed from the 
photometric redshift distribution by comparing the N(z) ob- 



tained in the center of the field with the equivalent one over a 
wider region of the same area, obtained under the same condi- 
tions from the z p hot point of view (same effective exposure time 
and number of filters), and used as a blank field. A real cluster 
or other structure should have appeared as an excess of galax- 
ies in the central region in comparison to the outer parts. In this 
exercise, we considered only objects with Ng > 1 i.e. objects 
that could not be excluded as galaxies without applying a cut 
in magnitude. Figure [6] displays the results found for the differ- 
ent fields. The histograms in this Figure display the difference 
between the redshift distribution within a ~ 140" radius region 
centered on the center of the image (N,„(z), black solid line), and 
the distribution within an outer ring, 140" < r < 200", (N OM ,(z), 
dashed black lines). Red solid lines show the positive difference 
between the two histograms, N,„(z) - N OM ,(z). Histograms were 
obtained with a Sz — 0.05 sampling step and smoothed with a 
5z = 0.15 sliding window. This window corresponded approxi- 
mately to the 1 <x uncertainty in the z P hot estimate for the faintest 
galaxies in the catalog. 

Where there was a distinct peak in N,„(z) - N ou/ (z) distribu- 
tion, we used this value to represent the "cluster redshift". We 
also computed 2D number density maps and cluster tomography 
(see § !5.3l below) to emphasize the reality of the clusters, in par- 
ticular for the uncertain cases. A summary of these results was 
provided in White et al. (2005). We note that the efficiency of the 
cluster-finding algorithm could be enhanced if the central region 
was centered on the cluster centroid instead of the center of the 
image. This ideal situation could be achieved in wider surveys. 

5.2. Spectroscopic versus photometric cluster redshifts 

Figure|6]summarizes the comparison between the z p hot and spec- 
troscopic redshifts for the different cluster fields. The photo- 
metric cluster redshift can be defined in different ways. Here 
we have adopted two different definitions, which are reported 
in Table [5] (1) the mean weighted value, computed from the 
excess peak, and (2) the redshift corresponding to the maxi- 
mum value in the N,„(z) - N OM/ (z) histogram. Table [5] provides a 
comparison between spectroscopic and photometric determina- 
tions of cluster redshifts for the low and high-z samples. Cluster 
redshifts in Table [5] and Fig. [6] correspond to the most promi- 
nent cluster identification when several clusters were present 
in the field (Mil vang- Jensen et al. 2008). Five fields were ex- 
cluded when computing the systematic deviation and dispersion: 
Cll 1 19-1 129 and C11238-1 144 in the low-z sample, because of 
incomplete photometry; and Cll 122-1 136, C11037-1243a and 
Cll 138-1 133 in the high-z sample, the first because of the lack 
of a clear cluster in the field, and the two others because their low 
cluster redshifts implied that B-band photometry was required to 
achieve an accurate z p hot (see also Milvang-Jensen et al. 2008). 
We note however the accurate photometric identification of the 
CI 1037- 1243a cluster, for which a poor determination was ex- 
pected. 

In general, the differences between photometric and spectro- 
scopic values are small, ranging from 5z ~0. 03-0. 04 at high- 
z to Sz ~0.05 at low-z. The dispersion is much lower than 
the approximate cut introduced by the spectroscopic preselec- 
tion {Zciuster ± 0.2). Therefore, there is no reason why a cluster 
should have been "missed" within the relevant redshift interval 
due to Zphot preselection. The systematic trend of lower quality 
in z p hot for the low-z sample was expected from the simulations 
presented in § 14.11 and also observed in the comparison with the 
spectroscopic sample in § 14.21 
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Figure 6. Photometric redshift distributions in the EDisCS fields, with the cluster redshift increasing from top to bottom and from 
left to right, for the low-z and high-z samples (first and second series respectively). The histograms display the following redshift 
distributions: N,„(z) (black thick solid lines), N oi( ,(z) (thin black lines), and the positive difference between N,„(z) - N OI „(z) (lowest 
histograms, thick red lines). A real cluster or other structure corresponds to a positive excess in the lower (red) histogram. Histograms 
were obtained with adz - 0.05 sampling step, smoothed by a 6z — 0.15 sliding window. Vertical solid and dashed lines indicate the 
spectroscopic and photometric values, respectively, adopted for the cluster redshift. All vertical scales are identical in the number 
of galaxies / Az = 0.05, from N(z)=0 to 100, apart from C11216-1201, for which the range is to 1 15. 



5.3. Cluster tomography 

We present results about cluster tomography along the line-of- 
sight in the different EDisCS fields. This highlights the capa- 
bility of Zphot in identifying and studying clusters of galaxies in 
deep photometric surveys. 

A local density estimate was derived at each point in the field 
using a grid with 5x = 6y =1". The density estimator was de- 
fined, according to the formalism introduced by Dressier (1980), 
to equal Z20 = 20/ (7rc/| ), where c/20 was the projected linear dis- 
tance to the 20-th closest neighbor. Photometric redshift slices 
of Az = 0.1 were used to cover the redshift range < z < 1 .4. 
Close neighbors were selected within z p hot ±0.1, centered on the 
redshift bin. The arbitrary choices of E„ and n=20 were justified 
by the typical richnesses of clusters and the z p hot accuracy. We 
tested other values of n ranging between 5 and 20 and achieved 
similar results for the cluster detection. Edge effects were cor- 
rected by using external fields in the main eight directions (i.e. 
X+, Y+, X-, Y-, X+Y+, X+Y-, X-Y+, X-Y-), containing differ- 
ent realizations of the same z p hot slice with xy coordinates ran- 



domly sorted. The density maps were smoothed with a Gaussian 
kernel of <x =5". 

Figures |7] and [8]display the projected number density maps 
obtained with this method, for the low and high-z samples, re- 
spectively. In these figures, clusters are presented with z c i uster in- 
creasing from top to bottom. Projected number densities are dis- 
played on a linear scale, for an arbitrary redshift step Az = 0.1, 
with isopleths corresponding to increasing number density bins, 
equally spaced with AE20= < ^20 >, starting at the mean £20 
within the redshift slice. The presence of a cluster along the line- 
of-sight is clearly seen in most cases, at least in all cases where a 
cluster was clearly detected in the field. We note that the < £20 > 
value is affected by the presence of a cluster in the correspond- 
ing redshift bin. The most significant contrast in the density map 
is found usually for the z c i uster ±0.1 redshift slice (rightmost col- 
umn of Fig. [7] and [8j. The position of the BCG (White et al. 
2005; Milvang-Jensen et al. 2008; Whiley et al. 2008; displayed 
by a blue cross in these Figures) coincides usually with the max- 
imum contrast in the density map for the most prominent clus- 
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Figure 7. Projected number density maps for the low-z sample, for different redshift slices with z c i us t er increasing from top to bottom. 
Projected number densities are displayed in a linear scale. Isopleths correspond to increasing number density bins with AX20 = 
< £20 >' starting at the mean X20 within the redshift slice. Thick frames highlight the redshift slices encompassing z c faster ± 0.1. 
The rightmost column displays the density map for the z c i uster + 0.1 redshift slice, where the maximum contrast in the density 
peak is usually reached. The position of the BCG for the most prominent clusters in the field are displayed by a blue cross. For 
C11059-1253, C11202-1224, and Cll 119-1 129, only two isopleths are displayed corresponding to 1 and 1.5 < E_o >■ Additional 
overdensities along the line of sight with detection levels exceeding 4<x are displayed by circles. High resolution figure available 
at http://www.ast.obs-mip.fr/users/roser/photoz_EDISCS/10644f8.pdf 



ters in a given field. The comparison between this maximum- 
contrast slice and the tomography with fixed arbitrary step (e.g. 
Az = 0.1, the step used in Fig. [7] and [8]) suggests that the opti- 
mal redshift step for detection of clusters along the line-of-sight 
should be close to the typical difference between photometric 
and spectroscopic z c i uster , in this case Az ~ 0.05. Structures sep- 
arated by less than Az ~ 0.05 in redshift space cannot be distin- 
guished by Zphot tomography (e.g. clusters Cll 138 (z=0.48) and 
C11138a (z=0.45) in the 1 138.2-1 133 field, and C11227 (z=0.63) 
and C11227a (z=0.58) in the 1227.9-1 138 field; Milvang-Jensen 
et al. 2008). 

For all clusters clearly identifiable in this sample (i.e. all 
clusters except C11059-1253, C11202-1224 and Cll 1 19-1 129), 
the maximum contrast is found to be at least ~ 3 - 4x < E20 > 
about z c i uster ±0. 1 . All clusters in our sample exhibit a significant 
> 3cr overdensity around z c i uster +Q. 1 , defined to be E20— < £20 > 
/cr(z), where cr(z) represents the standard deviation in the pro- 



jected number density within the redshift bin. For all apart from 
the three aforementioned clusters, the detection level exceeds 
4cr, and ranges between 6 and 9cr for C11216-1201, Cll 227- 
1138, C11411-1148, C11420-1236, C11040-1155, C11054-1245, 
C11138-1133,C11018-1211,and C11054-1146. No other signif- 
icant overdensities are found along the line-of-sight of density 
peaks exceeding 3x < Iqo >> which are typical values for rich 
clusters in this sample. However, several overdensities are found 
with slightly smaller values, around ~ 2 - 3x < E20 >> an d de- 
tection levels exceeding 4<x. These structures are identified by 
circles in Figures |7]and [8](see below). Their reliability is diffi- 
cult to assess with the presently available data. 

Several fields deserve further comment. More information 
concerning the spectroscopic identification of clusters in these 
fields can be found in the reference papers by Halliday et al. 
(2004) and Milvang-Jensen et al. (2008): 
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Figure 8. Projected number density maps for the high-z sample, for different redshift slices. Same com- 
ments as in Fig. [7] There is no clear BCG identified in Cll 122-1 136. High resolution figure available at 
http://www.ast.obs-mip.fr/users/roser/photoz_EDISCS/10644f9.pdf 



- C11301-1 139: The two clusters identified in this field with 
z(spectroscopic)=0.397 (Cll 301a) and 0.482 (C11301) are 
consistent with the two different > 3 < £20 > peaks observed 
by our tomography. 

- C11037-1243: The most distinct ~ 3 < S 2 o > over- 
density in this field corresponds to the "a" component at 
z(spectroscopic)=0.425, whereas the first identification was 
given at z(spectroscopic)=0.578 (z c ; BSter ±0.1 in Fig. [8]). Both 
structures are seen by tomography. 

- C11354-1230: the two clusters C11354 and C11354a 
(z =0.762 and z =0.595 respectively) are consistent with 
overdensities at the > 2 < Z20 > level. 

-Cll 103-1245: The prominent cluster detected in this field at 
z = 0.96 is clearly visible (> 8cr level) both in the density 
map and the z p hot distribution (Fig. [6]). The two components 
Cll 103a (z =0.626) and Cll 103b (z =0.703) can hardly be 
separated by tomography. 

- Other secondary peaks at z~ 0.9 - 1.1 are seen in the 
N(z p hot ) distribution of CI 1 1 8- 1 2 1 1 , Cll 30 1 - 1 1 39, CI 1420- 
1236, C11054-1245, and C11040-1155 (Fig. @. All can be 
associated with spatial overdensities of 2 - 3 < £20 > at the 



4 - 5cr level, apart from in the field of C11040-1 155, where 
no significant overdensity is found. 
- Several additional overdensities/structures at z~ 0.8 - 1.0 
are found in C1141 1-1 148, Cll 1 19-1 129, C11238-1 144, and 
C11216-1201, with 2-3 < Z20 > and detection levels ranging 
between 4 and 6 cr. Given the limited photometric coverage 
of Cll 1 19-1 129 and C11238-1 144, the detected overdensities 
in these fields are rather dubious. 

6. Cluster membership criteria 

The most unambiguous way to determine cluster membership is 
by means of accurate spectroscopic redshifts. Unfortunately, it is 
far too time-consuming to obtain high spectroscopic complete- 
ness in cluster member observations, even to relatively bright 
limits of / < 22 - 23. For this reason, it is necessary to de- 
velop membership criteria that rely solely on photometric data. 
To achieve many EDisCS science goals, such as study of lumi- 
nosity functions and cluster substructure, any method should: 1) 
retain >90% of cluster members, 2) reject an optimal number of 
non-members, and 3) measure a probability that a given galaxy is 
a cluster member. The first two criteria should be implemented 
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Table 5. Comparison between spectroscopic (z e /) and photomet- 
ric determinations of cluster redshifts in the low and high-z sam- 
ples. 



Cluster (low-z) 






r ^ j c\ ^phot 




^c/ ^ phot 






(1) 


(i) 


(2) 


(2) 


C11018-1211 


0.473 


0.575 


-0.102 


0.525 


-0.052 


C11059-1253 


0.456 


0.465 


-0.009 


0.478 


-0.022 


C11119-1129 1 


0.550 


0.544 


0.006 


0.475 


0.075 


CI 1202- 1224 


0.424 


0.556 


-0.132 


0.568 


-0.144 


C11232-1250 


0.541 


0.614 


-0.073 


0.575 


-0.034 


C11238-1144 1 


0.460 


0.548 


-0.088 


0.521 


-0.061 


C11301-1139 


0.482 


0.525 


-0.043 


0.524 


-0.042 


C11353-1137 


0.588 


0.579 


0.009 


0.573 


0.015 


C11411-1148 


0.520 


0.520 


0.000 


0.568 


-0.048 


C11420-1236 


0.496 


0.570 


-0.074 


0.569 


-0.073 


<A Z > 






-0.053 




-0.050 








± 0.050 




± 0.046 


<|AJ> 






0.055 




0.053 








± 0.048 




± 0.041 


Median A- 






-0.058 




-0.045 


Cluster (high-z) 




Zphol 


Zd . z ^ 


Z hot 


Zrf -Zphol 






(V 


c ^p o, 


(2) 


' (2) 


C11037-1243a' 


0.425 


0.461 


-0.036 


0.424 


0.001 


C11040-1155 


0.704 


0.635 


0.069 


0.624 


0.080 


C11054-1146 


0.697 


0.658 


0.039 


0.673 


0.024 


CI 1054- 1245 


0.750 


0.697 


0.053 


0.727 


0.023 


CI 11 03- 1245b 


0.703 


0.685 


0.018 


0.725 


-0.022 


Cll 122-1 136' 


0.640 


0.748 


-0.108 


0.773 


-0.133 


C11138-1133' 


0.479 


0.686 


-0.207 


0.720 


-0.241 


C11216-1201 


0.794 


0.747 


0.047 


0.725 


0.069 


C11227-1138 


0.635 


0.664 


-0.029 


0.625 


0.010 


C11354-1230 


0.762 


0.759 


0.003 


0.724 


0.038 


<A Z > 






0.028 




0.031 








± 0.033 




± 0.035 


<|A Z |> 






0.037 




0.038 








± 0.022 




± 0.026 


Median A- 






0.039 




0.024 



Notes - Two estimates for the photometric determination are provided: 
(1) mean weighted value, and (2) the redshift corresponding to the max- 
imum of the N(z) distribution. Clusters labeled with 1 were excluded 
when computing the systematic deviation and dispersion (see text). 

so that there is little dependence on the galaxy color, e.g. for 
Butcher-Oemler-type studies. While traditional methods of sta- 
tistical subtraction using "field" surveys of comparable depth of- 
fer a viable method to satisfy the first two criteria, they do not 
satisfy the third. For this reason, we developed an alternative 
method for membership determination based on our photomet- 
ric redshifts estimates. 

We present below how we use the photometric redshift prob- 
ability distribution P(z) to reject non-members from each cluster 
field. We describe the method that we developed and its cali- 
bration based on EDisCS spectroscopic redshifts. We discuss 
how this method can be extended to the entire magnitude-limited 
sample for a given cluster and outline its limitations. 

6. 1 . The method 

Traditionally, z^c-based methods for determining cluster mem- 
bership were based on a simple cut in redshift, such that a galaxy 
was considered to be a member if \z p hot - z c iust\ < ^Zthresh- One 
disadvantage of this method is that Az t hresh can be as high as 0.3 
(e.g. Toft et al. 2004), causing considerable field contamination 
to enter into the cluster sample. An additional disadvantage of 



the method is that it uses only the best-fit redshift in determining 
membership and ignores the information contained in the full 
redshift probability distribution P(z). 

Brunner & Lubin (2000) suggested an improved technique 
that used P{z) in determining galaxy membership. They assumed 
a Gaussian P{z) of width calculated from the comparison with 
Zspec , and defined the quantity 

fZdua+Az/2 

Pciust = P{z)dz, (1) 

Jz ri „ s ,-Az/2 

where Az was an interval defined around z c u<si that reflected 
the dispersion in the z p hot versus z spe c diagram, and z c iust was the 
spectroscopic redshift of the cluster. In Section I6.2l we explain in 
detail the method we used to calibrate a threshold value for P c i ust , 
P thresh, below which a galaxy would be considered to be a non- 
member. However, in reality, P{z) can be highly non-Gaussian, 
with multiple maxima and extended tails at large distances from 
the most likely solution. For this reason, it may not be optimal to 
assume the Gaussian approximation. We therefore extended the 
Brunner & Lubin method to use the full P(z) dataset calculated 
directly from the two z p h t codes. We tested the accuracy of our 
P(z) by comparing the confidence intervals derived from P(z) 
with the disagreement between z spe c and z P hot (see Figures [4] to 
[17] in §1?). The z sp ec fell within the 68% confidence limits on 
z P hot for ~ 68% of the galaxies and many of the galaxies with 
large \z p hot-z sp ec\ also had correspondingly large 68% confidence 
intervals. This gave us assurance that P(z) was accurate enough 
for our purposes. 

6.2. Calibrating from the spectroscopic sample 

We adopt the large and uniform EDisCS spectroscopic sample 
(Halliday et al. 2004; Milvang-Jensen et al. 2008) to calibrate 
the Pthresh that we use to reject non-members. 

We show in Figures [9] and [10] the P c i ust versus ZphorZrfust 
for all galaxies with secure z spe c measurements in our clusters, 
for the low and high redshift samples, respectively. Results pre- 
sented here were obtained for the GR code, but they are similar 
for Hyperz. As seen in Section |4] the z p hot accuracy as well as 
the fraction of catastrophic identifications in a given field depend 
on the spectral type of galaxies, although the difference between 
early and late types is smaller for the faintest galaxies in our sam- 
ple. Therefore, we have studied the reliability of the membership 
criteria for both early and late SED types using a cut in the rest- 
frame color which splits the sample into two equal halves of red 
(early) and blue (late) type galaxies. The color cuts are found to 
be (B - V) rest =0.79 and 0.67 for the low-z and high-z samples 
respectively. The left and right-hand panels in Figures [9] and [TO] 
provide results for objects of different spectral types. 

We note that the ratio of members to non-members increases 
as a function of P c iust- There are also few galaxies with z p hot 
~Z sp ec and very low P c iust values. This implies that there are not 
many members that would be rejected because their faint magni- 
tudes correspond to broad P(z) and lead to their rejection even if 
z P hot ~Z sp ec- It is important to consider, however, that the spectro- 
scopic sample consists of the brightest galaxies of probably the 
tightest P{z) values and that this behavior might not be similar at 
fainter magnitudes (see § |6.31 >. 

For the GR code, Figures QT] and [12] demonstrate how the 
retained fraction of members and rejected fraction of non- 
members depends on the P c i ust threshold P thresh for the low-z and 
high-z samples, respectively. In both cases, it is possible to de- 
fine a P^resh value such that >90% of confirmed cluster members 
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Figure 9. The integrated probability of being at z c iust ±0.1 vs. 
ZphorZciust for the low-z clusters. All galaxies with secure spec- 
troscopic redshifts were included. The different panels represent 
the blue and red halves of the sample in (B - V) rest . The typically 
lower values of blue galaxies is due to their broader probability 
distributions. The solid vertical line indicates z p hot -Zciust and the 
dashed lines indicate +0.1 in redshift. Horizontal error bars cor- 
respond to 68% confidence intervals. 



Table 6. Rejection thresholds 



Code 


Low-z 






High-z 






BVI 


BV1K 


BVIJK 


VRI 


VRIJK 


GR code 


0.350 


0.475 


0.300 


0.200 


0.050 


Hyperz 


0.150 


0.425 


0.425 


0.050 


0.400 



are retained with little dependence on rest-frame (or observed) 
color. 

Since the performance is similar using Hyperz, we use 
both codes jointly to provide the most efficient rejection. In 
Figs.[T8]and[TU we plot P cto (GR) versus P cto (Hyperz). While 
there is a large scatter, there is a definite correlation between 
the two probabilities, such that the majority of objects with 
low Pdust for one code also have a low P c i ust with the other 
code. A Spearman's Rank Correlation test on the distribution 
of PdusAGR) vs. P c /„ s/ (Hyperz) shows that there is higher than 
99.9% probability that these two variables are correlated. The 
same result is found for the low and high-z samples, when us- 
ing the full magnitude-limited samples or other subsamples re- 
stricted to the brightest galaxies. This implies that a joint re- 
jection is feasible. After extensive tests, we decided to reject 
galaxies if P c iust <Pthresh for either code; in these tests, the P thresh 
values for each code were determined separately, such that the 
highest rejection, independent of rest-frame color, was possible, 
while retaining >90% of the confirmed members. The adopted 
thresholds are summarized in Table [6] and the performance of 
these thresholds is summarized in Table Q 



Figure 10. Same as Figure|9]but for the high-z clusters. 



0.6 



0.4 



0.2 







• 4. ^ i; : 


.1 1 1 1 1 1 1 1 1 1 1 

••"mm. 




* A • 

6 


- ■ (B-V) re5t >0.79 

- ' (B-V) resl <0.79 
_ • all 

□ 

□ □ 
□ □ 


'a!. ■. 6 D 

DO* • 

Q 8 

□ □ □ 

□ □ ~© — 
□ o 2 






□ 




- 0°0 " 




" o 

A 









265 non-members _ 
178 members 





0.2 0.4 0.6 0.8 



thresh 

Figure 11. The fraction of retained members and rejected non- 
members as a function of the P t hresh an d as a function of (B - 
V)rest- This plot was created using all galaxies with secure Zspec 
measurements in clusters with B VIK photometry, clusters. Open 
symbols represent the fractions of rejected non-members, while 
solid points are the fractions of retained members. The color and 
shape of the points indicate the (B - V) res , cut applied. The solid 
horizontal line at 0.9 is included to guide the eye. This figure 
only shows the results for the GR code, as an illustration of the 
technique. It is not directly comparable to the numbers quoted in 
Tables [6] and [7] which utilize the combination of both the GR 
code and Hyperz. 
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Figure 12. Same as Figure[TT]but for clusters with VRIJK pho- 
tometry. 



Table 7. Retained and Rejected fraction in Spectroscopic 
Sample 





Low-z 

BVI 


BVIK 


BVIJK 


High-z 
VRI 


VRIJK 


ftnemb 
J retain 


0.95 


0.89 


0.90 


0.90 


0.89 


snonmemb 
J re ject 


0.53 


0.73 


0.88 


0.44 


0.53 



Because P(z) is broader for galaxies without NIR data, P c iust 
is also systematically lower and the P thresh, determined for galax- 
ies with NIR data, is no longer be applicable. To calibrate P thresh 
for galaxies without NIR data, we re-derived z P hot for the en- 
tire spectroscopic sample, excluding the NIR filters. We recal- 
ibrate P thresh an d summarize the performance and adopted cuts 
in Table [6] and [7] It is important to note that the performance of 
the rejection is different in areas with and without NIR data and 
the retained member population differs in the two regions. For 
this reason, we limit all studies using the photometric redshifts 
to those areas with NIR data. 

We checked how the effectiveness of the adopted P thresh var- 
ied across the sample. Because of the limited numbers of spec- 
troscopically observed objects per cluster, this was not possible 
on a cluster-by-cluster basis. Instead we split each of the high-z 
and low-z samples into two subsamples each and examined how 
the accepted and rejected fractions differed. The retained frac- 
tions of members ranges from 87-98% and the rejected fractions 
of non-members ranges from 50-60%. 



Figure 13. The apparent magnitude histograms for the spectro- 
scopic sample (solid blue histogram) versus the total sample of 
the same limiting magnitude (open black). For both samples at 
low-z (a and b) and high-z (c and d), two of the clusters with the 
most complete spectroscopic coverage are displayed. We note 
that the spectroscopic sample, even for clusters with the most 
extensive spectroscopy, is not complete at any magnitude limit 
and is biased towards brighter magnitudes. The magnitudes cor- 
respond to the /-band AUTO magnitudes from SExtractor. 



6.3. Pdust Threshold Performance in Magnitude Limited 
Samples 

We examine how applicable our adopted P t hresh> calibrated us- 
ing the spectroscopic subsample, is to the full magnitude-limited 
sample. In Fig. [13] we show the apparent magnitude distribution 
of the spectroscopic sample and the total photometric sample 
with the same magnitude limit, for two clusters in each redshift 
range with the widest spectroscopic coverage. It is clear from 
these plots that we are not spectroscopically complete at any 
magnitude and that the spectroscopic sample, as expected, is bi- 
ased towards brighter magnitudes. 

The spectroscopic target lists were not only constructed with 
a magnitude limit in mind, but also with an eye towards reducing 
the number of galaxies that had low probabilities of being at the 
cluster redshift, in addition to including some galaxies that were 
not formally present about Zciust- We examine in general how this 
preselection causes the P c iust distributions of the spectroscopic 
sample to differ from those in a magnitude-limited photomet- 
ric sample. In Fig. [20] we compare the histogram of the P c i ust 
values for the spectroscopic samples with those for the entire 
sample down to the same magnitude limits. The spectroscopic 
preselection manifests itself as an excess of high P c iust values 
and a deficit of low P c i u $t for the spectroscopic sample with re- 
spect to the photometric sample. Within the precision of our nu- 
merical routine, a KS-test gives 0% probability that these two 
distributions are drawn from the same parent distribution. This 
inherent bias implies that a certain P thresh removes a higher frac- 
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tion of galaxies in the photometric sample than was the case in 
the spectroscopic subsample. 

We examine how the spectroscopically calibrated P c i ust re- 
jection operates when applied to a magnitude-limited sample. 
We illustrate these results using the GR code, but the conclusions 
would be equivalent using Hyperz. In Figs.Q^JEUE] and|22] we 
plot the Pchst versus z p hot for two high-z and two low-z clusters. 
Each figure has panels that show how this distribution changes 
with apparent magnitude. As we move to fainter magnitude lim- 
its, many galaxies appear at all P c iust values. Those at high P c i ust 
do indeed fall close to Zdusu as predicted by the spectroscopic 
studies. Encouragingly, the galaxies with low P c i us t values fall 
systematically away from Zdust- In fact, even at the faintest mag- 
nitudes, there are very few galaxies within Zdust ±0- 1 that have 
Pciust < 0.2. We recall that this must not be the case: galaxies 
with Zphot "Zciust but broad P(z) distributions will have low prob- 
abilities of being at the cluster redshift, even though their best 
value lies around Zdust- This self-consistency implies that P c i ust 
is, in fact, providing us with a real indication of whether these 
faint objects are at the cluster redshift. This effect would be dif- 
ficult to reproduce by systematic errors in the photometric red- 
shifts because it appears for all clusters, regardless of their red- 
shift range or presence of NIR data. Nonetheless, there are some 
faint galaxies with z p hot ~ Zciust that may be rejected because of 
a broad P(z). Also, the tests presented in Sect. [43] demonstrate 
that the photometric redshift accuracy is expected to be lower 
for fainter galaxies, implying that there will be galaxies who are 
truly at Zduster but are scattered away from the cluster redshift. 
In Rudnick et al. 2009, we discuss how these effects may differ 
for red and blue galaxies and we present the implications for the 
study of the cluster galaxy luminosity function. 

The total fraction of galaxies rejected for each cluster as a 
function of / magnitude are presented in the columns 3-6 of 
Table[8] When our rejection criteria is applied, we reject 55-82% 
of the galaxies at / < 22 and 75-93% at / < 24.5. 




^phot 

Figure 14. A plot of P c Ust vs. z p i w t for C11018-1211 that com- 
pares the spectroscopic sample with the magnitude-limited sam- 
ple. Each tile corresponds to a different magnitude limit for the 
photometric sample. The full photometric sample is indicated by 
black dots. The spectroscopically confirmed cluster members are 
indicated by solid red triangles and the non-members by open 
blue squares. The solid vertical line indicates z p hot -Zdust and the 
dashed lines indicate the intervals +0.1 in redshift. We note that, 
at fainter magnitudes, the galaxies with low P c i ust values do not 
lie at Zphot ~ Zdust m large numbers, but rather are at different 
redshifts. 



Table 8. Total fraction of galaxies rejected as cluster members 
as a function of the / magnitude (columns 3 to 6 for / > 22, 23, 
24 and 24.5 respectively). 



Cluster Zci hi fa fit Ms 



C11018.5- 


1211 


0.4734 


0.68 


0.77 


0.85 


0.88 


C11037.5- 


1243a 


0.4252 


0.50 


0.58 


0.71 


0.77 


Cl 1040.4- 


1156 


0.7043 


0.68 


0.73 


0.83 


0.86 


C11054.2- 


1146 


0.6972 


0.56 


0.67 


0.78 


0.83 


C11054.4- 


1245 


0.7498 


0.58 


0.60 


0.69 


0.75 


C11059.1- 


1253 


0.4564 


0.68 


0.76 


0.84 


0.87 


C11103.4- 


1245b 


0.7031 


0.72 


0.77 


0.83 


0.86 


Cll 119.2- 


1129 


0.5500 


0.69 


0.75 


0.81 


0.83 


Cll 122.5- 


1136 


0.6397 


0.61 


0.66 


0.79 


0.84 


C11138.1- 


1133 


0.4796 


0.53 


0.62 


0.74 


0.79 


Cl 1202.4- 


1224 


0.4240 


0.80 


0.86 


0.91 


0.93 


C11216.4- 


1201 


0.7943 


0.62 


0.66 


0.74 


0.78 


C11227.5- 


1138 


0.6357 


0.58 


0.64 


0.73 


0.77 


C11232.3- 


1250 


0.5414 


0.67 


0.73 


0.82 


0.85 


C11238.3- 


1144 


0.4602 


0.65 


0.73 


0.81 


0.83 


C11301.4- 


1139 


0.4828 


0.71 


0.79 


0.87 


0.90 


C11353.0- 


1137 


0.5882 


0.77 


0.84 


0.89 


0.91 


C11354.1- 


1231 


0.7620 


0.64 


0.66 


0.75 


0.78 


C1141 LO- 


1148 


0.5195 


0.81 


0.84 


0.90 


0.92 


CI 1420.2- 


1236 


0.4962 


0.82 


0.87 


0.91 


0.93 



Notes - Only galaxies with NIR data were used in this calculation, 
except for Cll 119-1 129 and C11238-1144. For these two clusters, all 
galaxies were used because no NIR data exists. 
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Figure 15. Same as Figure [TJexcept for C11216-1201. 
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7. Discussion and Conclusions 

We have used two independent codes to compute photomet- 
ric redshifts: Hyperz and GR code. In general, the two codes 
yield rather similar results, either on a cluster-by-cluster basis 
or as a function of the filter set and spectral type, of typically 
cr(Az/(l + z)) ~ 0.05 to 0.06. Hyperz results are found to be 
slightly more accurate than GR's ones in general, by < 20% in 
c(Az/(l + z)). The quality achieved by both codes is consistent 
with the expectations derived from "ideal" simulations. An in- 
teresting trend is that the quality of both codes is highly corre- 
lated, in the sense that the highest and lowest quality results, in 
terms of <x(Az/(l + z)), and systematics are found for the same 
clusters. This trend cannot be due to the use of an incomplete 
or imperfect template set, as suggested by other authors (Ilbert 
et al. 2006), because in such a case, we should expect the same 
systematic behavior in all fields, given a filter set, as discussed 
in § 14.11 In contrast, different systematics are observed in the 
different fields, which are found to be almost equal for the two 
independent z P hot codes. This behavior suggests that the origin of 
the systematic errors is more likely to be associated with small 
residuals in the input photometry rather than the z p hot templates 
and codes. Indeed, small zero-point shifts of < 0.05 magnitudes 
cannot be excluded, in particular for the near-IR data. 

Photometric redshifts are found to be particularly useful in 
the identification and study of galaxy clusters in large surveys. 
The determination of cluster redshifts in the EDisCS fields us- 
ing a simple algorithm based on z p hot is highly accurate. Indeed, 
the differences between photometric and spectroscopic values 
are found to be small, typically ranging between 6z ~0. 03-0. 04 
in the high-z sample and 6z ~0.05 in the low-z sample. This 
is at least a factor ~ (1 + z) more accurate than the determina- 
tion of Zphot for individual galaxies. The accuracy is more sen- 
sitive to the filter set used rather than the redshift of the clus- 
ter. The systematic lower quality results for the low-z sample 
was somewhat expected from the simulations presented in § 14.11 
Tomography based on z p hot could be used in searches for clus- 
ters along the line-of-sight, using redshift steps optimized to be 
close in value to the typical difference between photometric and 
spectroscopic z c i uster to maximize the contrast between members 
and non-member galaxies (in this case, Az ~ 0.05). 

The cluster membership criterion presented in Sect. [6] has 
been used to extend the spectroscopic studies of cluster galaxies 
to fainter limits in magnitude (e.g. De Lucia et al. 2004, White 
et al. 2005, Clowe et al. 2006, Poggianti et al. 2006, De Lucia et 
al. 2007, Desai et al. 2007, Rudnick et al. 2009). 

In conclusion, photometric redshifts are useful tools for 
studying galaxy clusters. They enable efficient and complete 
pre-selection of cluster members for spectroscopy, allow accu- 
rate determinations of the cluster redshifts based on photometry 
alone, provide a means of determining cluster membership, es- 
pecially for bright sources, and can be used to search for galaxy 
clusters. 
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zfsDectroscoDic) — z(Dhotometric) 

Figure 16. Comparison between spectroscopic and photometric 
redshifts for the low-z sample, obtained with Hyperz (top panel) 
and GR (central panel) codes. Solid (red) circles, open (blue) 
circles, and crosses corresnond to obiects with pood rivne 1 V 



z(sT>ectroscoT>ic) — z(T)hotometric) 



Figure 17. (Comparison between spectroscopic and photomet- 
ric redshifts for the high-z sample, obtained with Hyperz (top 
panel) and GR (central panel) codes. Solid (red) circles, open 
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Table 9. z p hot accuracy derived from simulations for galaxies in the spectroscopic sample, based on the GR code. 



Clusters 


Redshift interval 


Galaxy type 


(A z > 


CTj. 


°~ z.MAD 


a(Az/(l+z)) 


1% 


g% 


Low-z 


0.3 < z plM < 1 


all 


0.008 


0.150 


0.104 


0.063 


0.9 


7.7 


High-z 


0.3 < Zpho, < 1 


all 


-0.072 


0.164 


0.167 


0.100 


10.6 


20.3 


High-z 


0.45 < Zpi, , < 1 


all 


-0.092 


0.158 


0.151 


0.088 


6.5 


21.5 


Low-z 


0.3 < Zpho, < 1 


E/S0 


-0.140 


0.129 


0.184 


0.110 


4.2 


6.3 




Sbc 


-0.077 


0.117 


0.136 


0.082 


0.8 


9.9 






Scd 


0.013 


0.076 


0.070 


0.041 


0.0 


4.3 






Im 


0.024 


0.077 


0.075 


0.045 


0.0 


5.9 






SB 


0.098 


0.117 


0.121 


0.076 


0.1 


9.8 


High-z 


0.3 < Zphot ^ 1 


E/S0 


-0.227 


0.153 


0.370 


0.226 


12.2 


31.6 




Sbc 


-0.101 


0.161 


0.232 


0.134 


17.3 


28.0 






Scd 


-0.025 


0.093 


0.086 


0.052 


2.7 


16.2 






Im 


0.024 


0.126 


0.123 


0.074 


5.9 


6.7 






SB 


-0.045 


0.063 


0.077 


0.047 


2.3 


13.5 



Notes - The accuracy reached for the different spectral types of galaxies is also presented. The information given is the same as in Tableland 
TablefT] for the current spectroscopic sample. 

Table 10. z p i,ot accuracy expected for the faintest galaxies in the EDisCS sample, based on the GR code. 



Clusters 


Redshift interval 


Galaxy type 


(A,) 




C z.MAD 


<r(Az/(l+z)) 


1% 


g% 


Low-z 


0.3 < Zpho, < 1 


all 


0.034 


0.136 


0.124 


0.082 


3.1 


0.0 


High-z 


0.3 < z^o, < 1 


all 


0.003 


0.131 


0.134 


0.081 


3.4 


2.9 


High-z 


0.45 < Zpho, < 1 


all 


0.016 


0.124 


0.130 


0.078 


3.4 


3.0 


Low-z 


0.3 < zpho, < 1 


E/S0 


0.015 


0.123 


0.118 


0.078 


1.8 


0.0 




Sbc 


0.043 


0.128 


0.124 


0.081 


0.7 


0.0 






Scd 


0.036 


0.145 


0.112 


0.075 


6.4 


0.0 






Im 


0.093 


0.179 


0.153 


0.103 


8.9 


0.0 






SB 


0.041 


0.145 


0.129 


0.083 


6.7 


0.0 


High-z 


0.3 < Zpho, < 1 


E/S0 


0.014 


0.129 


0.149 


0.088 


1.2 


1.6 






Sbc 


-0.010 


0.125 


0.128 


0.080 


0.9 


1.3 






Scd 


-0.008 


0.134 


0.114 


0.067 


4.0 


4.7 






Im 


0.012 


0.124 


0.122 


0.073 


5.4 


0.0 






SB 


0.011 


0.140 


0.134 


0.092 


13.6 


11.7 



Notes - The accuracy reached for the different spectral types of galaxies is also presented. The information given is the same as in Tableland 
Table[2] for the current spectroscopic sample. 
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Figure 18. A comparison of P c iust computed by the two different photometric redshift codes for the clusters with B VIK photometry. 
The black open circles represent all galaxies with / < 24.5. The solid red squares are spectroscopically confirmed members and the 
solid green triangles are spectroscopically confirmed non-members. There is a broad correlation of P c iust between the two codes, 
such that the overwhelming number of galaxies have low P c imt f° r both methods with a smooth distribution of higher P c iust galaxies 
extending to the upper right, coincident with the confirmed members. The blue lines indicate the P t hresh values for each code that 
were determined to reject jointly the largest number of non-members, while retaining at least 90% of the confirmed members. 
Objects to the left of the vertical blue line or beneath the horizontal blue line are flagged as interlopers. 
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Figure 19. Same as Figure[l8]but for the clusters with VRJJK photometry. 
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Figure 20. The histogram of P c i ust values for the photometric sample (solid line) versus that of the spectroscopic sample (dotted 
line). Panel (a) is for the high-z sample at / < 22.5 and panel (b) is for the low-z sample at / < 22. In each case, the histograms 
represent the fraction of the total number of objects over all P c iust- 
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Figure 21. Same as Figure fl4lexcept for C11420-1236. 
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Figure 22. Same as Figure [T4l except for C11040-1 155. 



